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We study the dependence of the electronic density of states (DOS) on the distance from the 
boundary for a ferromagnet/superconductor bilayer. We calculate the electron density of states in 
such structure taking into account the two-band model of the ferromagnet (FM) with conducting s 
and localized d electrons and a simple s-wave superconductor (SC). It is demonstrated that due to 
the electron s-d scattering in the ferromagnetic layer in the third order of s-d scattering parameter 
the oscillation of the density of states has larger period and more drastic decrease in comparison 
with the oscillation period for the electron density of states in the zero order. 

PACS numbers: 74.45.-l-c, 74.78.Fk 

As it well known superconductivity creates a gap in the electronic density of states (DOS) near the Fermi energy Ep. 
In the case of ferromagnet/superconductor multistructures at the SC/FM boundary the Cooper pair wave function 
extends from the SC into the FM layer. The internal exchange field in a ferromagnet results in a strong suppression of 
the superconducting order parameter. Damped oscillatory dependence of the Cooper pair wave function in FM hints 
that a similar damped oscillatory behavior may be expected for the variation of the DOS due to proximity effect. For 
SC/FM bilayer, for example, this question has been studied in It was shown that the DOS oscillations in FM 
near the boundary are described by the expression 

Nf{e = 0) « NiO) (^1 - lexpi^^)cosi^)^ , (1) 

where is the characteristic length of the superconducting correlations decay in FM layer, x is the distance from the 
boundary. 

In most theoretical papers P, the nonmonotonic behavior of the superconducting order parameter was studied 
in the so-called dirty limit and for low energic^ that is using the quasiclassical Usadel equations in the context of 
one-band model of ferromagnetic metal. In ref.yf the more general Eilenberger equations for a larger energy range are 
discussed. However, due to some simplifications, in both cases there are evident limitations which could be avoided 
by solving the system of full Gor'kov equations for the normal and anomalous Green functions. Moreover, when 
calculating the DOS in 3d ferromagnetic metal one should take into account the s-d electron scattering which is the 
main scattering mechanism responsible for the kinetic properties of 3d metal due to the large DOS of d electrons at 
the Fermi surface. 

As it was substantiated above we study the influence of the s-d scattering of the conducting electrons on the DOS 
in the ferromagnetic layer taking into account the two-band model of the FM in contact with the SC using Green 
function approach. The FM is assumed to be a 3d ferromagnetic metal with two types of electrons - conducting s 
electrons and almost localized d electrons. The SC is a simple s-wavc superconductor. 

We consider the planar semi-infinite geometry with an FM to the left of z=0 and an SC to the right. The Hamiltonian 
in the F layer is 

+ EaE„/rfr5(p-p„)7«[V'+^V'^+V'rC], ^ ' 

in the S-layer is 

H^J2 [ dr^'t'^sr^ + [ dr[A(r)V.+>J+ + c.c], (3) 

where ipi^'''^^ ^nd ipcr''^'^ are the creation and annihilation operators for s(d) electrons with spin cr, A(r) is the 
superconducting order parameter, ^s(d) are the dispersion functions for the s(d) electrons, 7„ is the electron 
random s-d scattering parameter. In the calculations we use the averaged over impurity configurations value 
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of 7 = 7" and (A7) = 7,^ — ^ . The normal and anomalous Green functions of the one-band model are 
G^'*(a;i, 2:2) = — (Tt^^ {xijij}^'^ {X2)) and F'^'^[xi,X2) = —(Tr^^'^{xi)ip^'^{x2)), where x = (r, r) is a four-component 
vector and the creation and annihilation operators are associated with s conducting electrons. The two-band model 
requires us to take into account similar Green functions of d electrons and the Green functions responsible for the s-d 
electron scattering. Then in the FM layer with s-d scattering the full Green function can be written as the following 
matrix in the s-d space: 
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in the SC layer the same matrix has the form: 
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In the spin space, each of the functions G"^ and F*-' is also the two by two matrix: 
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With the accoimt of the term responsible for the s-d electron scattering in (1) the system of the Gor'kov equations for 
normal and anomalous Green functions become rather complicated. Such a system is no longer linear and requires 
the perturbation solution over the parameter 7 which is small in comparison with the electron kinetic energy. 

In the zero order on 7 we have the anomalous function F only of the s electrons due to the boundary conditions 
with the S-layer. In order to obtain all other functions in (3) the next approximations have to be investigated. In the 
second order on the s-d electron scattering parameter the Gor'kov equations contain F"^^ function which is connected 
with F^* by means of 7. Anomalous F'""^ and F'^'^ functions in the F-laycr originate from their connection with G'^'^ in 
higher order on 7. We are interested in the third order on the s-d scattering parameter 7 as in this case the solution 
of the Gor'kov equations has the form of plane waves with the exponents responsible for the change of the electron 
density of states oscillation period due to s-d scattering. 



We consider the full Green function as the sum of the components of the corresponding order on 7: G = G 



(2) 



andF = F(o)+F(i)+F, 



(0) 



^(2) . .- By carrying out the Fourier transformation in the plane perpendicular 
to the z axes (using the quasi momentum representation in the xy plane) in the third order on 7 we get the following 
system of Gor'kov equations in the F-layer: 



(z, z' , w) — 7 ■ Gj.| {z, z', uj) = 0, 
7 ■ G\f'{z, z', a;) + . G^f ^(z, z', u;) ^ T^., 
-El-Flf^\z,z',u;)+j-Flf'%z,z',u;) = 0, 
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where kp^^^^ is the Fermi impulse of s(d) electrons with the spin up (down); kpsup the electron Fermi impulse in the 
S-layer; mg(^j^ the mass of s(d) electrons; ta the energy parameter; ep the Fermi energy; xq the impurity concentration. 
The functions G^^\G^'^\f'^°\ and F'^) are considered to be known from the previous order on 7. In particular the 
normal diagonal (z — z') Green function of the zero order G*-*^^ has the form 

G(°)"' (z, z) = consti + const2e^'''^^'>^^ , where (9) 

the coefficients 

consti, 2=consti, 2(fc(o)i, /c(o)2, fcs, t^, A, {/^^Y) 
are the functions of 

/c(o)i - ^^4" -X^ + 2»m,a. - 2m,xo(A7)2Gff (z, z), (10) 
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superconducting order parameter A, energy parameter w, and s-d electron scattering parameter 7. One can see that 
in the zero order on 7 the electron density of states of s electrons near the boundary with the S-layer has a rapid 
oscillation with the period proportional to T^/k'^p^ and exponential decaying to the bulk value at large z. In this 
case, the conducting electrons undergo the ordinary and Andreev reflection at the border F/S which leads to this 
rapid oscillations. Another type of oscillations due to s-d electron scattering appears in higher order on 7 when the 
relationship between the normal s electron Green function and anomalous takes an explicit form. Then we get the 
superposition of the ordinary rapid oscillations due to Andreev reflection on the interface and slower ones with the 
period proportional to the reversed difference of the Fermi momenta of s electrons with spin up and spin down. 
In the S-layer the system of Gorkov equations for the normal and anomalous Green functions keeps its usual form: 
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The Green functions G'^'^ and F'^'^ in S-layers are identically zeros. 

In the third order of 7 the diagonal function G^3^(z,z,x) which serves to calculate the DOS in the FM layer 
represents the following sum of the plane waves with the different wave vectors: 
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where the coefflcients are the functions of the bi = ^(Aj)'^ fi{ki, k2, k^, k^, kl 
As one can see the first term in (7) contains the exponent with the argument proportional to the doubled difference 
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of the Fermi impulses of s electrons with the opposite spin directions {kp 
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That is, this term gives the 



increase of the characteristic period of oscillations due to s-d scattering. The effective exponential decaying is defined 
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/(2TOsa;o(A7) /mG||^j^-)) is the free path of the s electron with the 
is the mass of the s electron, xo is the impurity density, G'^^in is the 



by the exponent e ^'tt+It)^ where l^^^^ = 
spin up(down) in ferromagnetic layer, ms 

Green function of d electrons. In order to calculate the DOS we use the imaginary part of the diagonal Green 
function: p''(z) = j^Im J xdxG'"' {z , z , x) ■ The graphical dependencies of the DOS of s electron on the distance from 
the boundary in the ferromagnetic layer in the zero order and in the third order for the different values of 7 and 
fixed (A7)2 are presented in Fig. 1(a) and 1(b), correspondingly. The oscillations of the third order term in the s 
electron density of states have larger period and more drastic decrease in comparison with the oscillation period for 
the electron density of states in the zero order. This fact evidences the influence of the s-d electron scattering of 
conducting electrons on their density of states in the ferromagnetic layer. 

On the FM side, the proximity effect induces the superconducting order parameter. Increasing the distance from the 
interface it displays a damped oscillation and changes sign from positive to negative. Usually the states corresponding 
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Figure 1: (a) DOS as a function of the distance z from the boundary in FM in the zero approximation, (b) The additional 
term to the DOS as a function of the distance z from the boundary in FM for = and three different values of s-d electron 
scattering parameter 7 (k\^ = 1.2 A"\ fe^^ = 0.42 (A7)^ = 0.1 eV^). 




Figure 2: Energy variation of the DOS in FM for different distances z from the boundary. Zero energy corresponds to the 
Fermi level, energy gap A = 1.4 meV, exchange field /i = 10 meV. 



to a positive sign of superconducting order parameter are called the " state" and those corresponding to a negative 
sign of order parameter the "tt state". In ref.Q Guoya et al. develop a quantum-statistical approach based on the 
McMillan and BTK theories to calculate the DOS dependence on the energy. They show that the DOS in FM displays 
a maximum at the energy-gap edge and a minimum at the Fermi level exhibiting the SC-like shape for the "0" state, 
and flipped shape for the "tt" state. 

Figure 2 shows the result of our calculations of the DOS in FM using the Gor'kov equations. As it may be seen the 
DOS exhibits the same as in ref. Q peak-dip behavior near the SC gap for "0" and "tt" states. With increasing the 
distance from the interface, the SC-like behavior of the DOS in F-layer gradually disappears. Such two different DOS 
shapes in FM near the FM/SC boundary for the "0" and "tt" states have been observed experimentally by Kontos et 
al. [a| in the measurements of the DOS by planar-tunnelling spectroscopy in Al/Al203/PdNi/Nb functions. 

On the SC side the superconducting order parameter is diminished near the FM/SC interface (see Fig. 3) and 
owing to the proximity effect the superconductivity in the SC near the interface also becomes gapless. Far from the 
boundary the DOS in SC takes it usual shape as in bulk superconductor. 
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Figure 3: Energy variation of the DOS in the SC near (solid line) and far (dotted line) from the boundary. 
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